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Abstract: We study the ground-state octet baryon masses and sigma terms using 
the covariant baryon chiral perturbation theory (ChPT) with the extended-on-mass-shell 
(EOMS) renormalization scheme up to next-to-next-to-next-to-leading order (N 3 LO). By 
adjusting the available 19 low-energy constants (LECs), a reasonable fit of the nj = 2+1 lat- 
tice quantum chromo dynamics (LQCD) results from the PACS-CS, LHPC, HSC, QCDSF- 
UKQCD and NPLQCD collaborations is achieved. Finite-volume corrections to the lattice 
data are calculated self-consistently. Our study shows that N 3 LO BChPT describes better 
the light quark mass evolution of the lattice data than the NNLO BChPT does and the 
various lattice simulations seem to be consistent with each other. We also predict the pion 
and strangeness sigma terms of the octet baryons using the LECs determined in the fit of 
their masses. The predicted pion- and strangeness-nucleon sigma terms are a n N = 43(1) (6) 
MeV and a sN = 126(24) (54) MeV, respectively. 
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1 Introduction 

Studies of baryon spectroscopy play an important role in understanding the nonpertur- 
bative nature of quantum chromodynamics (QCD). In the past few years, Lattice QCD 
simulations [1, 2] has made remarkable progress in studies of non-perturbative strong- 
interaction physics. Recently, the lowest-lying baryon spectrum, composed of up, down 
and strange quarks, has been studied by various LQCD collaborations [3-11] . However, 
because these calculations adopt different lattice setup and all of them lead to the same 
continuum theory, it is crucial to test whether the results for baryon masses are consistent 
with each other [12]. On the other hand, since lattice QCD simulations are performed in 
a finite hypercube and with larger than physical light quarks masses [13], the final results 
can only be obtained by extrapolating to the physical point (chiral extrapolation) and in- 
finite space-time (finite volume corrections). Chiral perturbation theory (ChPT) provides 
a useful framework to perform such extrapolations and to study the induced uncertainties. 
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ChPT, as the low energy effective field theory of QCD, is based on effective Lagrangian 
techniques with an expansion in powers of external momenta and light quark masses, con- 
strained by chiral symmetry and its breaking pattern [14-25]. Its applications in the 
mesonic sector have been rather successful at least in the two flavor sector of u and d 
quarks. But the extension to the one-baryon sector turns out to be non-trivial. Because 
baryon masses do not vanish in the chiral limit, a systematic power-counting (PC) is ab- 
sent [17]. In order to restore the chiral power-counting, the so-called Heavy-Baryon (HB) 
ChPT was first proposed by Jenkins and Manohar [26] , considering baryons as heavy static 
sources. Although this approach provides a strict power-counting, the heavy baryon ex- 
pansion is non-relativistic, which might lead to pathologies in certain cases and is found 
to converge rather slowly in the three flavor sector of u, d, and s quarks. Later, covariant 
BChPT implementing a consistent PC with different renormalization methods have been 
developed, such as the infrared (IR) [27] and the extended-on-mass-shell (EOMS) [28, 29] 
renormalization schemes. In addition to the afore-mentioned dimensional renormalization 
schemes (MS and its derivatives), to speed up the convergence of BChPT, other renor- 
malization/regularization schemes are also proposed, e.g., the cutoff scheme [30], the finite 
range regulator (FRR) method [31-33], and the partial summation approach [34]. 

In the past decades, the ground-state (g.s.) octet baryon masses have been studied 
extensively [33-48]. It is found that SU(3) HBChPT converges rather slowly [49]. Further- 
more, most calculations are performed only up to NNLO because of the many unknown 
low-energy constants (LECs) at N 3 LO except those of Refs. [38, 39, 41, 46-48]. Regarding 
chiral extrapolations, Young and Thomas [33] obtained very good results using the FRR 
scheme up to NNLO by fitting the LHPC [7] and PACS-CS [5] lattice data. In Ref. [49], 
we applied the NNLO EOMS-BChPT to analyze the same lattice data and found that the 
EOMS-BChPT can provide a better description of lattice data and is more suitable for 
chiral extrapolation purposes than HBChPT and NLO BChPT. Recently, using a partial 
summation scheme up to N 3 LO, Semke and Lutz [25, 46, 48] found that the BMW [4], 
HSC [8], PACS-CS [5], LHPC [7], and QCDSF-UKQCD [10] lattice results can be well 
described. 

On the other hand, up to now, a simultaneous description of all the nj = 2 + 1 lattice 
data with finite volume effects taken into account self-consistently is still missing. 1 Such 
a study is necessary for a clarification of the convergence problem and for testing the 
consistency between different lattice simulations. Furthermore, it also provides a good 
opportunity to determine/constrain the many unknown LECs of BChPT at N 3 LO. 

In this work we study the g.s. octet baryon masses and sigma terms using the EOMS- 
BChPT up to N 3 LO. Finite-volume corrections to the lattice data are calculated self- 
consistently and are found to be important in order that a good fit of the lattice data can 
be achieved. Unlike Refs. [33, 46], the virtual decuplet contributions will not be explicitly 
included, because its effects can not be disentangled from those of virtual octet baryons 
due to the large number of unknown LECs. 2 In order to fix all the 19 LECs and test 

1 In Ref. [48], Semeke and Lutz showed that their partial summation approach can reproduce the results 
of the HSC and QCDSF-UKQCD collaborations by fitting the BMW, PACS and LHP data. 

2 However, their inclusion might be important for the determination of baryon sigma terms, as having 
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the consistency of current lattice calculations, we perform a simultaneous fit of all the 
publicly available n f = 2 + 1 LQCD data from the PACS-CS [5], LHPC [7], HSC [8], 
QCDSF-UKQCD [10] and NPLQCD [11] collaborations. 

The paper is organized as follows. In Sec. 2 we collect the relevant terms of the chiral 
effective Lagrangians up to N 3 LO, which involve 19 unknown LECs. The explicit results 
of the g.s. octet baryon masses are provided in Sec. 3, where the EOMS renormalization 
scheme is also briefly explained. Sec. 4 focuses on the determination of the unknown LECs 
by fitting all the LQCD data and the chiral extrapolation of the octet baryon masses. We 
calculate the pion- and strangeness-baryon sigma terms in Sec. 5. A brief summary is given 
in Sec. 6. 



2 Power counting and effective Lagrangians 

2.1 Power counting 

In ChPT, the power counting (PC) provides a systematic organization of the effective La- 
grangians and the corresponding loop diagrams within a perturbative expansion in powers 
of (p I A x sB) nchPT ) where p is a small momentum or scale and AchPT the chiral symmetry 
breaking scale. In the one-baryon sector, the chiral order, n^px, °f a diagram with L 
loops is calculated as 

n ChPT = 4L - 2N^ - N B + ^ kV k , (2.1) 

k 

where N^(Nb) is the number of internal meson (baryon) propagators, and Vk is the num- 
ber of vertices from /cth-order Lagrangians. However, as mentioned in the introduction, 
in the covariant BChPT this systematic power counting is lost. That is to say, in the 
calculation of a loop diagram one may find analytical terms with a chiral order lower than 
that determined by Eq. (2.1)[16]. In order to recover the power counting, we adopt the so 
called extended-on-mass-shell renormalization scheme. In this heme, the lower-order power 
counting breaking pieces of the loop results are systematically absorbed into the available 
counter-terms [27, 29] . A detailed discussion of the EOMS renormalization scheme will be 
presented in Sec. 3.2. 

2.2 Chiral Lagrangians 

In this subsection, we collect the relevant chiral Lagrangians for the calculation of octet 
baryon masses in the three-flavor sector of u, d and s quarks up to N 3 LO. The Lagrangians 
can be written as the sum of a mesonic part and a meson-baryon part: 

L cS - + L (j} + L 4>B + L (j}B + L 4>B + L^g, {2.2) 

where the subscript (i) denotes the corresponding chiral order 0(q l ), (j) = (ir, K, rj) 
represent the pseudoscalar Nambu-Goldstone boson fields, B = (N, A, S, H) the g.s. 
octet baryons. 



been stressed in Ref. [50] and will be the subject of a forthcoming work. 
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2.2.1 Meson Lagrangians 

The lowest-order meson Lagrangian is given by 

4^ = ^l(D,U(D,U^) + 3 (x[ /t + C7x t) ! (2 .3) 

where Fa, is the pseudoscalar decay constant in the chiral limit, (X) stands for the trace 
in flavor space, x = ^BqM accounts for explicit chiral symmetry breaking with Bq = 
— (0\qq\0)/F^, M. = diag(m;, mi, m s ), where we assumed perfect isospin symmetry, m u = 
rrid = mi- The 3x3 unimodular, unitary matrix U collects the pseudoscalar fields 

U{4>) = u 2 (<p) = exp (if) , (2.4) 



with 



a=l 



f^ + fa vr+ K + \ 



(2.5) 



V K~ K° -fa J 

Under SU(3) L xSU(3) K , U(x) transforms as U -> U' = LUR\ with L,R G SU(3) L ,«. 
The most general meson Lagrangian at 0{q i ) has the following form [16]: 

4 4) = L x [{D li U{D^)f^L 2 {D li U{D v U)^){D^U{D v U)^) 

+L 3 (D^U(D»UyD u D(D u Uy) + U{D^U(D^){ X U ] + U X ] ) 
+L 5 {D lM U(D^U)\ X U ] + U X ] )) + L 6 [( X tf + U X ] )f 
+L 7 [( X rf - U X ] )] 2 + L s {Ux ] U X ] + XU ] XU ] ) 
-iL 9 (f*D»U(D»Uy + f^D^D'U) + L 1Q (Uf 

+Hi{f*fg' + fau /l" ) + ^2(xX f ), (2-6) 



A" 



where /£" = - ^V" - i[r**, r"] and /£" = W - W - Z"] with = v M + a 
= v fj, — a n with and the external vector and axial currents. 
The LECs Li 's are scale-dependent and absorb the infinities generated by the one-loop 
graphs. 

2.2.2 Meson-Baryon Lagrangians 

The effective meson-baryon Lagrangians contains terms of even and odd chiral orders, 

reft _ , r (?) , ^(3) , A*) / 9 7 n 

<j>B (f>B ' 4>B ' <f>B <f>B' 

The lowest-order meson-baryon Lagrangian is 

= - mo)5> + ^(BT^K^Ji), (2.8) 

where mo denotes the baryon mass in the chiral limit, and the constants D and F are 
the axial-vector coupling constants, which are determined from the baryon semi-leptonic 
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decays. The traceless 3x3 matrix B contains the lowest-lying octet baryon fields 



B a ^a 



^ y/2 

(2=1 



V 



(2.9) 



Under SU(3)lxSU(3)^, B transforms as any matter field, B -> B' = KBK\ with 
K(U, L, R) the compensator field representing an element of the conserved subgroup SU(3)y • 
In equation (2.8), the covariant derivative of the baryon field is defined 



D/J.B = dfj,B + [r^, B], 



and Up the axial current defined as 



(2.10) 
(2.11) 

(2.12) 



where u = yU. 

The meson-baryon Lagrangian at order 0(q 2 ) can be written as 

,(2)' 



c 



(2) 



r (2, sb) rV 



(2.13) 



This separation is motivated by the fact that the first part appears in the tree and loop 
graphs, whereas the latter only contributes via loops. The explicit chiral symmetry breaking 
part reads as 

) = bo(x + )(BB) + b D/F (B[ X+ ,B] ± ), (2.14) 

where bo, bo, and bp are LECs, and x+ = u^x u ^ + u\^u. For the latter part, we take the 
same form as in Ref. [51]: 

C% = hiBlu^iu^B]]) +6 2 (%,K,5}}) 
+b z {B{u^[u^B]}) +b A {BB){u^) 
+ib 5 ((5K, [u v ^D v B\\) - (B%[u v , K, 7 ^]]) 

+i& 6 ((fl[<K,7 M A,J3}]> - (bKK,K, 7/1 B]}) 

+i& 7 ((B{/,{/, 7(1 D„B}}) - (flKK,K,7,J3}})) 

+i6 8 ([B^D V B) - (B*D vlll B)) {u»u u ) + ■■■ , (2.15) 

where &i ... 4 have dimension mass -1 and 65,... ,8 have dimension mass -2 . If one works for a 
set of fixed quark masses (e.g., Ref. [42]), all terms with one or two covariant derivatives 
can be absorbed in the structures proportional to 6i,...,4. However, for our purposes, we 
need to retain all the terms because they lead to different quark mass dependencies. 
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The third chiral order Lagrangian does not contribute to the baryon masses and the 
fourth-order effective Lagrangian relevant to our calculations is [38]: 

4fl = d 1 (B[ X+ ,[x+,B}]) + d 2 (B[ X+ ,{ X+ ,B}]) 
+d 3 (B{ X+ , {X+,B}}) + d A {B X+ ){ X+ B) 
+d 5 (B[ X+ ,B])( X+ )+d 7 (BB)( X+ ) 2 

+d 8 (BB)( X l). (2.16) 
In total, up to N 3 LO, we have 19 LECs: mo, bo, b B , bp, 6i_g and c?i_5 j 7 , g. 

3 Octet baryon masses 

In this section, we evaluate the octet baryon masses up to fourth order using the covariant 
BChPT, supplemented by the extended-on-mass-shell (EOMS) renormalization scheme. 

3.1 Baryon self-energy up to N 3 LO 

The physical baryon mass is defined at the baryon pole, ft = m B , in the two-point function 
of the baryon field ip B (x) 

S (x) = -i(0\T[Mx)MVW) = 

ft -mo - E{ft) 

where mo is the baryon mass in the chiral limit and S($) corresponds to the baryon self- 
energy 

m B - m - E(ft = m B ) = 0, => m B = m + E(ft = m B ). (3.2) 

The leading contribution to the self-energy, S a = Tn§ , is of order 0(q 2 ) [Fig. 1(a)]. The 
self-energy Y,^ = m^ of the one-loop diagram [Fig. 1(b)] is of order 0(g 3 ). One tree 
diagram contribution from Ct B [Fig. 1(c)] and two loop diagrams [Figs. l(d,e)] are of order 

0(q 4 ), m™ = S c + Erf + E e . We remark that due to parity conservation, there are no first 
order contributions. The baryon mass up to fourth order in the chiral expansion can be 
expressed 

m B = mo + m^+mf+mf. (3.3) 

The tree diagrams Fig. l(a,c) can be calculated straightforwardly. The corresponding 
results are shown below. The three one-loop diagrams Figs. 1(b), 1(d) and 1(e) yield, 
generically, 

Gb = 1 I Wf h5 ft-t-m + ie h5 k*-M* + ie> (3 ' 4) 

/d 4 k 1 
(2^ & W ^7-} fc2 _ M| + , V (3-5) 

(2 ) f d 4 k „ ( 1 \ 2 „ 1 



G e = i m> B > / -—j ^75 - — j, — ^7570 , f2 , ■ ■ ( 3 - 6 ) 

(27r) 4 \p — ft — mo + iej k — Mi + ie 
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Figure 1. Fcynman diagrams contributing to the octet-baryon masses up to 0(<? 4 ) in the EOMS- 
BChPT. The solid lines correspond to octet-baryons and dashed lines refer to Goldstone bosons. 
The black boxes (diamonds) indicate second (fourth) order couplings. The solid dot (circle-cross) 
indicates an insertion from the dimension one (two) meson-baryon Lagrangians. Wave function 
renormalization diagrams are not explicitly shown but included in the calculation. 

where Mm represents the mass of a Nambu-Goldstone boson. Here, we want to mention 
that the integral of Eq. (3.4) has been calculated in Ref. [44], and the results can be 
found there. The above loop functions contain power-counting breaking (PCB) terms 
and therefore additional steps need to be taken to conserve a proper chiral power-counting 
scheme. Among the different approaches, the EOMS scheme has been shown to be superior 
to heavy-baryon or infrared approaches for a number of observables (at least at the one- 
loop level) [44, 52-54]. 

3.2 The EOMS renormalization scheme 

The essence of the EOMS scheme is to perform additional subtraction of PCB pieces 
beyond the MS or MS renormalization scheme. That is possible because the PCB terms 
appearing in a loop calculation are analytical and can always be removed by redefining the 
corresponding LECs. This is equivalent to removing the finite PCB pieces directly from 
the loop results. For an explicit calculation, the PC can be restored in two slightly different 
ways: (1) one can first perform the loop calculation analytically, and then remove the PCB 
terms, or (2) one can first perform an expansion in terms of the inverse heavy-baryon mass, 
1/rriB, calculate the PCB terms, and then subtract them from the full results. It should 
be noted that the first approach may not always work because one may not be able to 
calculate the loop diagrams analytically, but since the PCB terms are finite and analytical, 
the second prescription should always work. 

In the present study, we have explicitly checked that all the PCB terms appearing in 
our loop calculation can be removed by redefining the LECs introduced in the previous 
subsection. 

3.3 The mass formulas 

After calculating all the Feynman diagrams shown in Fig. 1 and subtracting the PCB terms 
using the EOMS renormalization scheme, we obtain the full expressions of the octet baryon 
masses up to N 3 LO. 
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At 0(q 2 ) the tree level contribution provides the leading-order (LO) SU(3)-breaking 
corrections to the chiral limit octet baryon mass 



• (2) " E 4>J, (3.7) 

=7T, K 



where the coefficients £g . are listed in Table 1 

Table 1. Coefficients of the NLO contribution to the self-energy of the octet baryons (Eq. (3.7)). 





N AS 




t(a) 

*(«) 

^B,K 


-(26 + 46 F ) ^(36 - 26 D ) -(26 + 46 D ) 
- (46 + 46 D - 46 F ) ^ (66 + 86^ ) - 46 


-(26 - 46 F ) 
-(460 + 4615+46^) 



At 0{q^) diagram Fig. 1(b) gives the NLO SU(3)-breaking corrections to the baryon 
masses 

=7r, K, r) 



>B=Jj^ry E (3.8) 



The coefficients are given in Table 2, and the corresponding loop functions can be 
found in Ref . [44] . 

Table 2. Coefficients of the NNLO contribution to the self-energy of the octet baryons (Eq. (3.8)). 





N 


A 


s 




Ab) 

Ab) 
>>B,K 
Ab) 


§(L> + F) 2 

l(5D 2 - 6DF + 9F 2 ) 
l(D-3Ff 


2D 2 

§(D 2 + 9F 2 ) 
ID 2 


l(D 2 + 6F 2 ) 
2{D 2 + F 2 ) 

id 2 

3^ 


l(D-F) 2 

i(5L> 2 + 6DF + 9F 2 ) 
\{D + W) 2 



The NNLO SU(3)-breaking corrections to the octet baryon masses are 

+74^ E [&$B$ x) m + ^>r ] m) + ^h£*\m,) 

^ r ^ ^ =7r , K, V 

+71^2 E eS^-^&w (3-9) 



(4tt 2 F^) 2 

B'=N, A, s, s 

The first three terms of Eq. (3.9) are the tree contributions of diagram Fig. 1(c), and the 
corresponding coefficients £5^., £b\k can ^ e f° un d i n Table 3. The next term is the 

contribution from the tadpole diagram Fig. 1(d) and the Clebsch-Gordan coefficients are 
listed in Table 4. The last term is from the one-loop diagram of Fig. 1(e), together with the 
wave function renormalization diagrams not shown, and Cbb' s can ^ e f° un< i in Table 5. 
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The loop functions, which we use the MS renormalization scheme to cancel the divergent 
pieces and the EOMS renormalization scheme to remove the PCB terms, are written as 

,2 



4 d)1) (M^ = Ml 
Hf 2) {M^) = Ml 

r (d,3) 



1 + ln 

1 +ln 

M 4 



M l 

1 + ln 



Ml 



(3.10) 
(3.11) 
(3.12) 



h b]b'( M 4>) 



IMl 



mhUml-Ml 



6m 2 , (m^ - m$) - M|(2m£ ; - m$) 

2^(2) _(2) 



,(2) _(2)> 



-Mi 



arccos ■ 



{2) (2) 3^- m ^)-M^- m ^) M 
3(m^ — m^, ) H n — z m 



rat. 



2niQ 

(2) m (2^ A/f 2 
2 



???, 



+ (mg ) +mg ) )ln^ 



m 



21 





(3.13) 



(2) (2) 

where the m B and m^/ are the corresponding LO SU(3) corrections to the octet baryon 
masses given in Eq. (3.7). 

Table 3. Coefficients of the N 3 LO tree contribution to the sclf-cncrgy of the octet baryons 
(Eq. (3.9)). 





N 


A 










-4(4di +2d 5 +d 7 + 3d 8 ) 


-4(4d 3 + |d 4 + d7 + 3d 8 ) 


-4(4d 3 


+ d 7 + 3d 6 


) -4(4di-2d 5 + d 7 + 3d 8 ) 


*(<0 


-16 (di - d 2 + d 3 
— d 5 + d-7 + d 8 ) 


-16(§d 3 + §d 4 + d7 + d 8 ) 


-I6(d 7 


+ d 8 ) 


-16 (di +d 2 +d 3 
+d 5 + d 7 + d a ) 


M 


8 (4di - 2d 2 - d 5 
-2d 7 + 2d 8 ) 


16(§d 3 + ^d 4 -d 7 + d 8 ) 


-16(d 7 


-de) 


8 (4di + 2d 2 + d 5 
-2d 7 + 2d 8 ) 



At N 3 LO, a replacement of the meson masses by their 0(q 4 ) counterparts in 



??? 



(2) 
B 



generates N 3 LO contributions to mS . The corresponding Nambu-Goldstone boson masses 
up to 0(q 4 ) can be found in Ref. [16], which read as 



„ 2 f Ml 2 (Ml 2 

Ml A = Ml A 1 + ^44. In' 



M 



V .2 



327r 2 F 2 



96tt 2 F 2 



In 



M 2 



16 



^ [(2m, + m s )B Q (2Ll - L\) + m,B (2LS " L s)] 



(3.14) 



M|, 4 



9 f M «2 







M 2 



16 



(2m z + m s )B (2Ll - L\) + -(m, + m s )£? (2^ - Lg) 



(3.15) 
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Table 4. Coefficients of the tadpole contribution to the self-energy of the octet baryons (Eq. (3.9)). 





AT 


A 


S 




Mi) 
Ml) 

<*B,K 


-3(26 + b D + b F )ml 
-2(46 + 3b D - b F )m% 
-i [8(6 + 6b> - fcF)mi 

d L * J / -TV 

-(26o + 3& D -56 F )m2] 


-2(36 + 6 D )m2 
-§(66 + 5b D )m 2 K 
-§ \4(3b + 4b D )m 2 K 
-(360 + 76^)77x2] 


-6(60 + 6^)7772 
-4(26 + b D )m 2 K 
-f Inborn 2 , 
+(6 D - 60)7772] 


-3(26 + 6 D - 6 F )m2 
-2(460 + 36^ + 6^)777^ 
-i |"8(6 + 6/j + 6 F )mi 
-(26o + 36c + 56 F )m2] 


Mfi) 
^B,n 

Ma) 

S>B,K 
M,2) 
SB,?? 


3(6i + 6 2 + b 3 + 26 4 ) 
2(36i + 36 2 - 6 3 + 46 4 ) 
|(96i + 62-363 + 664) 


2(26 2 + 36 4 ) 
|(96! + 6 2 + 664) 
2(26 2 + 64) 


2(46! + 26 2 + 36 4 ) 
4(6i + 6 2 + 26 4 ) 
|(26 2 + 36 4 ) 


3(6i + 6 2 - 63 + 26 4 ) 
2(36j+ 36 2 + 63 + 464) 
i(96i + 6 2 + 36 3 + 66 4 ) 


M,3) 

Sb,ot 

M3) 

^B,K 

MS) 


6(65 + 6 6 + 6 7 + 26 8 ) 
4(36 5 - 6 6 + 36 7 + 46 8 ) 
|(96 5 - 36 6 + 6 7 + 66 8 ) 


4(26 7 + 36 8 ) 
|(96, 5 + 6 7 + 66 8 ) 
4(26 7 + 6 8 ) 


4(46, 5 + 26 7 + 36 8 ) 
8(65 + 6 7 + 26 8 ) 
|(26 7 + 36 8 ) 


6(65 - 6 6 + 6 7 + 26 8 ) 
4(36,5 + 6 6 + 36 7 + 46 8 ) 
|(96 5 + 36 6 + 6 7 + 66 8 ) 


Table 5. Coefficients of the loop contributions (Fig. lc) to the self-energy of the octet baryons 
(Eq. (3.9)). 


N 


A 


S 






= 41.= 

£nn v — Tui® ~ 3F) 2 Caai; = 
£nak = T2~( D + 3F ) 2 41 = 

<r( e ) 3/n FX 2 — 


-\{D + ZF) 2 41 = 
\& 41 = 

d 2 41 = 
i^-3^) 2 41 = 

M 

SSHA' — 


= h(D-F) 2 41 

1 n 2 C ( e ) 
3^ SHSA' 

2F 2 4 e i 

ln2 <r„„ 

3"^ S==r7 - 

±{D + F) 2 


= ^(D-3F) 2 
= \{D + F) 2 
= \{D-F) 2 
= ^D + 3F) 2 
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and 



The empirical values of 2Lg — L\ 
3Lj + Lg = —0.15 x 10 -3 are taken from the latest global fit [55], which are evaluated at 
the renormalization scale = 0.77 GeV. 3 To be consistent with our renormalization scale 



3 It should be noted that the uncertainties of the L\ are quite large. Because the effects of their 
contributions are found to be small, we do not take into account the uncertainties of these LECs in our fit 



of the LQCD mass data. 
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used for the one-baryon sector, we have re-evaluated the L£'s at \i = 1 GeV. 
4 Chiral extrapolation of octet baryon masses 

In this section, we study light quark mass dependence of the g.s. octet baryon masses using 
the N 3 LO EOMS-BChPT mass formulas by fitting the n f = 2 + 1 lattice QCD simulation 
results. As mentioned above, there are 19 unknown LECs, which cannot be fully determined 
by the lattice data of a single LQCD simulation. Therefore, we decide to fit all the publicly 
available lattice results of the g.s. octet baryon masses obtained by different collaborations, 
including the PACS-CS [5], the LHPC [7], the HSC [8], the QCDSF-UKQCD [10] and the 
NLPQCD [11] collaborations. In performing such a study, we can test the consistency 
of all these lattice simulations, which used totally different setup. Furthermore, because 
different calculations are performed with different lattice spacing, the uncertainties of our 
results also incorporate (partly) discretization effects. 

4.1 2+1 flavor LQCD data of octet baryon masses 

As has been stated in the introduction, the lattice calculations are performed with larger 
than physical light quark masses and finite volumes. In Fig. 2 , we show the lattice simula- 
tion points of the PACS-CS, LHPC, HSC, QCDSF-UKQCD and NPLQCD collaborations 
in the (2M^- — M%) -M% plane and in the L—M% plane. There are only one point for the 
NPLQCD simulation at m n = 389 MeV on the left panel, because the main purpose of the 
NPLQCD collaboration is to study the finite-volume effects on the baryon masses. The 
large range of light pion masses provides an opportunity to explore the applicability of 
ChPT for extrapolation of baryon masses. Although the light u/d quark masses adopted 
are always larger than their physical counterpart, the strange quark masses vary from col- 
laboration to collaboration: those of the PACS-CS and LHPC collaborations are larger 
than the physical one; those of the HSC and NPLQCD groups are a bit smaller, whereas 
those of the QCDSF-UKQCD collaboration are all lighter than the physical one. 

In the L—M% plane, it is seen that the PACS-CS and LHPC simulations adopt a single 
value of lattice volume; the HSC and QCDSF-UKQCD simulations use two different lattice 
volumes and the NPLQCD simulations are performed with four different lattice volumes in 
order to study the finite- volume effects on the octet baryon masses. Many of the simulations 
are still performed with M^L from 3 to 5 and with M& larger than 300 MeV. As a result, 
finite- volume corrections may not be negligible (see, e.g., Ref. [56]). In our study, we will 
take into account finite- volume corrections in a self-consist way as in Refs. [44, 56]. It 
should be stressed that taking into account finite- volume corrections is not only needed to 
obtain a good fit of the LQCD data but also can constrain better the values of some of the 
LECs, since they contribute to finite-volume corrections. 

Except for varying light- and strange-quark masses and lattice size, lattice simulations 
can adopt different fermion actions, all of which are believed to lead to the same continuum 
theory. Therefore, it's crucial to test whether all these simulation results are consistent 
with each other [12]. 
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Figure 2. (Color online). Landscape of LQCD simulations of the ground-state octet baryon masses: 
the PACS-CS (red circles), LHPC (blue squares), QCDSF-UKQCD (green diamonds), HSC (yellow 
upper triangles) and NPLQCD (sky-blue pentagons) Collaborations in the 2Mj C — M% vs. M% plane 
(left panel) and in the L vs. M% plane (right panel). The star denotes the physical point with the 
physical light- and strange-quark masses. 



In Appendix A, we tabulate the octet baryon masses of the PACS-CS, LHPC, HSC, 
QCDSF-UKQCD and NPLQCD collaborations. The numbers are given in physical units 
using either the lattice scale specified in the original publications [5, 7, 8, 11] or the method 
of ratios such as QCDSF-UKQCD [10]. N 3 LO BChPT is not supposed to be able to describe 
all the lattice simulations, which can have very large light quark/strange quark masses and 
quite small volumes. Recent studies show that m, < 0.4 — 0.5 GeV may be already at the 
limit where N 3 LO BChPT can work. Therefore, to reduce contributions from higher order 
terms in the chiral expansion, we choose the lattice simulations with M 2 < 0.25 GeV 2 . To 
minimize finite volume corrections, we further require the lattice data to have M^L > 4 . 
There are only 11 lattice data sets satisfying both requirements: three from PACS-CS, two 
from LHPC, only one from QCDSF-UKQCD, two from HSC and three from NPLQCD. We 
denote these lattice sets with stars in Tables 9 - 13 of Appendix A. For comparison, we also 
choose a large data set with M 2 < 0.5 GeV 2 , M% < 0.7 GeV 2 and M^L > 3. There are 
26 lattice data sets left and they are listed in Table 9 (six sets from PACS-CS), Table 10 
(five sets from LHPC), Table 12 (eight sets from QCDSF-UKQCD), Table 11 (four sets 
from HSC) and Table 13 (four sets from NPLQCD). In the following, we denote the above 
11 lattice sets as Set-I, and the 26 sets as Set-II. It should be mentioned that we are well 
aware that N 3 LO BChPT may not be able to describe the larger data set, which is indeed 
the case as we will see in the following subsection. 

4.2 Chiral extrapolation of octet baryon masses 

We proceed to fit the LQCD results of octet baryon masses by using the N 3 LO EOMS- 
BChPT mass formulas Eq. (3.3). For the meson decay constant, we use = 0.0871 
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GeV [57]. 4 For the baryon axial coupling constants D and F, we use the SU(6) relation 
F = 2/3D. Together with D + F = 1.26 as determined from nuclear beta decay, one then 
has D = 0.8 and F = 0.46. 5 The renormalization scale [i is set at 1 GeV, as in Ref. [44]. 

Set-I and Set-II contain 11 and 26 sets of data from the PACS-CS, LHPC, HSC, 
QCDSF-UKQCD and NPLQCD collaborations, respectively. Apparently, the results of 
different collaborations are not correlated with each other. But the data form the same 
collaboration are partially correlated by the uncertainties propagated from the determina- 
tion of the lattice spacing. Therefore, in order to correctly calculate the x 2 we incorporate 
the inverse of the resulting correlation matrix Cij = OiOjbij + AajAaj for each lattice en- 
semble (see Ref. [44]) , where the crj are the lattice statistical errors and the Aaj are the 
fully-correlated errors propagated from the determination of aj. The finite- volume correc- 
tions to the baryon masses are consistently calculated in the EOMS-BChPT framework, 
and the details will be reported in Ref. [58]. 

We perform a \ 2 fit to the lattice data and the physical octet baryon masses by varying 
the 19 LECs. The so-obtained values of the LECs from the best fits are listed in Table 6. 
We remark that in the following, the experimental octet baryon masses are always included 
in the fit unless otherwise stated. This way, the prediction of the sigma terms of the octet 
baryons and the fitted values of LECs are better constrained. 

For the sake of comparison, we have fitted Set-I using the NLO 6 and NNLO EOMS- 
BChPT. The values of the LECs bo, bo, bp, and ttiq are tabulated in Table 6. An order- 
by-order improvement is clearly seen, with decreasing % 2 /d.o.f. with each increasing order. 
Apparently, only using the 0(q 3 ) chiral expansion, we cannot describe simultaneously the 
lattice data from the five collaborations. The corresponding x 2 /d.o.f . is about 8.6. On the 
other hand, in the N 3 LO fit of lattice data Set-I and experimental octet baryon masses, 
X 2 /d.o.f. = 1.0. In addition, the values of the fitted LECs (named Fit I) all look very 
natural. Especially, the baryon mass in the chiral limit m,Q = 880 MeV seems to be 
consistent with the SU(2)-BChPT value [59, 60]. The fit to data Set-II yields a x 2 / d -°- f - 
about 1.6 and the fitted LECs look similar to those from Fit I except 62, b§, 67 and b$. 7 
The increased % 2 /d.o.f. indicates that data set II is a bit beyond the application region of 
N 3 LO BChPT. 

In Fig. 3, setting the strange-quark mass to its physical value, we plot the light-quark 
mass evolution of N, A, X and E as functions of M % using the LECs from Table 6. We can 
see that the NNLO fitting results are more curved and do not describe well lattice data 

4 In principle at N 3 LO one can use either the chiral limit value, F$ — 0.0871 GeV, or the SU(3) averaged 
value, F4, = 1.17-FV with = 0.0924 GeV as in Ref. [44]. The difference is of higher chiral order. In 
practice, we found that at N 3 LO the results are not sensitive to these two options while at NNLO the 
SU(3) averaged value is more preferred by the LQCD data. 

5 We have allowed D and F to vary in the fits and found that the optimal values determined by the 
Lattice data are consistent with the phenomenological values. 

6 Since at this order ChPT does not generate any finite-volume corrections, we have subtracted from the 
LQCD data the finite- volume corrections calculated by N 3 LO EOMS-BChPT with the LECs determined 
from the corresponding best fit. 

7 It should be noted that some of the LECs are correlated and other data except the masses are needed 
to fully determine the LECs. 
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Table 6. Values of the LECs from the best fit to the LQCD data and the experimental data at 
0(q 2 ), 0(q 3 ), and 0(q 4 ). The estimator for the fits x 2 /d.o.f. is also given (see text for details). 







Set-I 




Set-II 


Fit - 0(q 2 ) 


Fit - 0(q 3 ) 


Fit I - 0{q 4 ) 


Fit II - 0{q A ) 


m [MeVl 


900(6) 


767(6) 


880(22) 


868(12) 


b fGeV" 1 ] 


-0.273(6) 


-0.886(5) 


-0.609(19) 


-0.714(21) 


b D [GeV" 1 ] 

L J 


0.0506(17) 


0.0482(17) 


0.225(34) 


0.222(20) 


b F [GeV" 1 ] 


-0.179(1) 


-0.514(1) 


-0.404(27) 


-0.428(12) 


61 [GeV" 1 ] 






0.550(44) 


0.515(132) 


62 [GeV" 1 ] 






-0.706(99) 


0.148(48) 


63 [GeV" 1 ] 






-0.674(115) 


-0.663(155) 


fe 4 [GeV" 1 ] 






-0.843(81) 


-0.868(105) 


fe 5 [GeV" 2 ] 






-0.555(144) 


-0.643(246) 


b 6 [GeV" 2 ] 






0.160(95) 


-0.268(334) 


fe 7 [GeV" 2 ] 






1.98(18) 


0.176(72) 


b$ [GeV" 2 ] 






0.473(65) 


-0.0694(1638) 


di [GeV" 3 ] 








0.0340(143) 


0.0345(134) 


d 2 [GeV" 3 ] 






0.296(53) 


0.374(21) 


d 3 [GeV" 3 ] 






0.0431(304) 


0.00499(1817) 


d 4 [GeV" 3 ] 






0.234(67) 


0.267(34) 


d 5 [GeV" 3 ] 






-0.328(60) 


-0.445(26) 


d 7 [GeV" 3 ] 






-0.0358(269) 


-0.183(12) 


d 8 [GeV" 3 ] 






-0.107(32) 


-0.307(21) 


X 2 /d.o.f. 


11.8 


8.6 


1.0 


1.6 



Set-I. On the contrary the two N 3 LO fits, named Fit I and Fit II, both can give a good 
description of lattice data Set I. The rather linear dependence of the lattice data on m 2 at 
large light quark masses, which are exhibited both by the lattice data [7] and reported by 
other groups, is clearly seen. 

Next, we would like to test the predicting power of N 3 LO EOMS-BChPT. To do that, 
we take the same fitting strategy as before but remove the experimental values from the 
fit. The fits of data Set-I and Set-II yield a x 2 /d.o.f. about 0.93 and 1.6, respectively. The 
extrapolated results at the physical point are listed in Table 7, in comparison with the 
experimental masses [12]. The error bars are the uncertainties propagated from the fitted 
LECs. The N 3 LO extrapolations are in better agreement with the measured masses than 
the NNLO extrapolations. Surprisingly, the extrapolated results using the LECs from fit 
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Figure 3. (Color online). The lowest-lying baryon octet masses as functions of the pion mass. The 
two bands correspond to the best 0(q A ) fit to lattice data Set-I and Set II, and the dot-dashed lines 
and dashed lines are the best NLO and NNLO fits to lattice data Set-I. In obtaining the ChPT 
results, the strangeness quark mass has been set to its physical value. The lattice numbers are 
projected ones with N 3 LO BChPT with the LECs determined from the best fit to set-I and their 
strange quark mass is also set to the physical value. 



Table 7. Extrapolated baryon octet masses to the physical point of the 0(q 2 ), 0(q 3 ), and 0(q A ) 
best fits. The error is only statistical. 



B 


Exp. [12] 


Fit I 


Fit II 


Fit-I (NLO) 


Fit-I (NNLO) 


N 


0.940(2) 


0.865(26) 


0.936(10) 


0.987(6) 


0.876(3) 


A 


1.116(1) 


1.076(18) 


1.129(8) 


1.118(7) 


1.067(4) 


S 


1.193(5) 


1.155(17) 


1.193(7) 


1.177(6) 


1.126(4) 




1.318(4) 


1.294(23) 


1.310(8) 


1.279(6) 


1.288(4) 



II seem to be closer to the data than those using the LECs from Fit I, though the latter 
fit yields a smaller x 2 /d.o.f.. The same is true for the NLO ChPT fit. It yields the largest 
X 2 /d.o.f. but the extrapolations are even closer to the physical values than the N 3 LO Set-I 
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Figure 4. (Color online). The PACS-CS, LHPC, QCDSF-UKQCD and HSC data in comparison 
to the BChPT fits as functions of the pion mass. The curves from bottom to top are for N, A, 
£ and S, respectively. The bands in each panel are the Fit-II results at the 68% confidence level 
with M\ = a + bM%, where a and b are determined by fitting M\ as a linear function of M% for 
the corresponding lattice ensemble (see Appendix B). It should be noted that because one of the 
PACS-CS data (denoted by the solid pink point) has a small strange-quark mass, it is not included 
in the M\ fit. Finite-volume corrections have been subtracted from the lattice data using the 
corresponding BChPT fit. 

fit. This fact shows that the LQCD mass data with ~ 0.3 — 0.5 GeV alone do not 
constrain well the LECs. 

As mentioned above, the Fit-II is a global fit to the 26 lattice data sets with a 
X 2 /d.o.f. = 1.6. In Fig. 4, we show its description of the lattice data from the PACS-CS, 
LHPC, QCDSF-UKQCD and HSC collaborations, respectively. 8 The baryon masses are 
plotted as functions of M%, with the kaon mass calculated using = a + bM% with a 
and b determined from the lattice data for each ensemble (see Appendix B for details), 
and the lattice data are all extrapolated to infinite space-time. It is clear that our fitting 
results can give a reasonable description of the lattice data, with % 2 /d.o.f. equals 0.93, 
1.38, 0.93, 0.92 for each collaboration's data, respectively. It should be noted that Fig. 4 
is only for demonstration purposes because we should not trust the fit-II results, which 

8 We do not show the NPLQCD data because they are obtained at a single pion mass. 
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might be misleading because of the inclusion of lattice data beyond the application region 
of N 3 LO BChPT. 



5 Pion- and strangeness-baryon sigma terms 

In this section, we evaluate the pion- and strangeness-sigma terms for the octet baryons at 
the physical point using the N 3 LO BChPT results of the previous section. Baryon sigma 
terms are quantities of great importance. In particular, the nucleon-sigma term is a vital 
quantity in understanding the composition of the nucleon mass and the strangeness content 
of the nucleon. An accurate knowledge of the strangeness-nucleon sigma term is also useful 
for dark matter searches [61]. ChPT, with its LECs fixed by the LQCD mass data, can 
make predictions for sigma terms [62-64]. 

The sigma terms are the scalar form factors of baryons at zero recoil. In this work, we 
calculate all the octet baryon sigma terms o+b, &sB for B = N, A, E, H , by use of the 
Feynman-Hellmann theorem, which states: 

8Mb 

a nB = m l {B(p)\uu + dd\B(p)) = m i-^- (5.1) 

dM B 

a sB = m s {B(p)\ss\B(p)) = m s — . (5.2) 

8m s 

where mi = (m u + m^)/2. Using leading-order ChPT, the quark masses can be expressed 
by the pseudoscalar masses, with mi = M%/(2B ) and m s = (2Mj c — M%)/(2B ). Other 
related quantities, which often appear in the literature, including the strangeness content 
(i/b) and the so-called "dimensionless sigma terms" (fiB, Isb), are also calculated: 

_ 2{B(p)\ss\B(p)) _ mi 2a sB 
(B(p)\uu + dd\B(p)) m s a nB 
mi{B(p)\uu + dd\B(p)) a nB 

= Mb = M~b (5 - 4) 

_ m s (B(p)\ss\B(p)) _ a sB , , 

fsB ~ Mb " M~ B ( } 

Using the previous Fit-I LECs, we obtain the pion- and strangeness sigma terms cr n B, 
a S B for all the octet baryons, and the corresponding strangeness content i/b, "dimensionless 
sigma terms" Jib, fsB- The results are shown in Table 8. 

In Fig. 5, we compare our results of the sigma terms with the nj = 2 + 1 1QCD results 
of Refs. [62, 65] and the latest results obtained in the NNLO HBChPT with the FRR 
scheme [63] and the N 3 LO covariant BChPT with the partial summation scheme [64]. Our 
calculated baryon sigma terms are consistent with those of others, except for a n ~. and ct s n'- 
Our a n ~, is a bit smaller and ct s n = 126(24) (54) is slightly larger than the other results,. 

The nucleon pion-sigma term at the physical point, a n N = 43(1) (6) MeV, is in rea- 
sonable agreement with the determination in the study of the old ir — N scattering data 
[66] of a n N = 45 ± 8 MeV, but smaller than the central value of the more recent study, 
= 59 ± 7 MeV [67]. Our a n N is also in agreement with the recent lattice result 
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Figure 5. (Color online) Pion and strangeness sigma terms of the octet baryons. The solid 
circles represent our results. The empty circles are the results from the BMW [62], QCDSF [65] 
collaborations and from the BChPT studies with the FRR scheme [63] and the partial summation 
scheme [64]. 



of the QCDSF collaboration (a nN = 38(12) MeV [68], cj^tv = 37(8)(6) MeV [69]) and 
BMW Collaboration (ct^n = 39 (4) ^ 8 MeV) [62] within uncertainties, and is consistent 
with the HBChPT result of Ref. [63]. On the other hand, our value is larger than the 
QCDSF result (ct^n = 31(3) (4) MeV) [65] and the N 3 LO BChPT results using the partial 
summation scheme (ovjv = 32(1) MeV) [64], but slightly smaller than the JLQCD result 
(a wN = 50(4.5) MeV) [70]. 

It should be mentioned that the central value of the nucleon sigma term a n N in our 



Table 8. The sigma-tcrms, the strangeness content and the "dimcnsionless sigma terms" of the 
octet baryons at the physical point. The first error is statistical and the second one is systematic, 
estimated by taking half the difference between the N 3 LO result and the NNLO result. 



a nB [MeV] a sB [MeV] 



VB 



fl 



IB 



fsB 



N 43(1)(6) 126(24)(54) 

A 19(1)(7) 269(23)(66) 

£ 18(2)(6) 296(21)(50) 

S 4(2) (3) 397(22) (56) 



0.244(47)(110) 
1.179(118)(522) 
1.369(180)(512) 
$.263(4157) (6306) 



0.0457(11)(64) 
0.0170(9) (63) 
0.0151(17)(50) 
0.00303(152) (228) 



0.134(26)(57) 
0.241(21)(59) 
0.248(18)(42) 
0.301(17)(42) 
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present calculation is smaller than that of Ref. [44], where the result is obtained with NNLO 
EOMS-BChPT by fitting the PACS-CS data. Inclusion of virtual decuplet contributions 
may have some non- negligible effects on the predicted baryon sigma terms. At NNLO, it 
is found that the inclusion of virtual decuplet baryons can increase the pion-nucleon sigma 
term while decrease the strangeness-nucleon sigma term. It is interesting to check whether 
such effects still exist at N 3 LO. 

6 Summary 

We have studied the lowest-lying octet baryon masses with the EOMS BChPT up to N 3 LO. 
The unknown low-energy constants are determined by a simultaneous fit of the latest 2 + 1 
flavor LQCD simulations from the PACS-CS, LHPC, HSC, QCDSF-UKQCD and NPLQCD 
collaborations. Finite- volume corrections are calculated self-consistently. It is shown that 
the eleven lattice data sets with M 2 < 0.25 GeV 2 and M^L > 4 (0 = tt, K, rj) can be 
fitted with a x 2 /d.o.f. = 1.0. Including more lattice data with larger pion masses or smaller 
volumes deteriorates the fit a bit but still yields a reasonable x 2 /d.o.f = 1.6. 

Our studies confirm that covariant BChPT in the three flavor sector converges as 
expected, i.e., relatively slowly as dictated by mR-/A x sB but with clear improvement order 
by order, at least concerning the octet baryon masses. A successful simultaneous fit of all 
the latest 2+1 flavor LQCD simulations indicates that the LQCD results are consistent 
with each other, though their setups are quite different. 

Applying the Feynman-Hellmann theorem, we have calculated the sigma terms of the 
octet baryons. In comparison with the results obtained in previous studies, we conclude 
that to obtain the sigma terms with an accuracy of a few percent, more works are still 
needed. 
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A Nf = 2 + 1 Lattice QCD simulation results 

We briefly summarize some key ingredients of the LQCD simulations of the PACS-CS [5], 
LHPC[7], HSC [8], QCDSF-UKQCD[10] and NPLQCD [11] collaborations, which are rele- 
vant to our study. In addition, we tabulate the simulated octet baryon masses in physical 
units, which satisfy m 2 < 0.5 GeV 2 , m 2 K < 0.7 GeV 2 , and M^L > 3. 

PACS-CS [5] 

The PACS-CS collaboration employs the nonperturbatively 0(a)-improved Wilson quark 
action and the Iwasaki gauge action. Numerical simulations are carried out at the lattice 
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spacing of a = 0.0907(14) fm, on a 32 3 x 64 lattice with the use of the domain-decomposed 
HMC algorithm to reduce the up-down quark mass, which is about 3 MeV. For the strange 
quark part they improve the PHMC algorithm with the UV-filtering procedure. Their 
simulation points cover from 701 MeV to 156 MeV, but the lightest point has a m^L ~ 2.9, 
which might induce large finite volume corrections. 



Table 9. Masses of the pseudoscalar mesons and the octet baryons (in units of MeV) obtained by 
the PACS-CS collaboration (TABLE III of Ref. [5].) The first number in the parentheses is the 
statistical uncertainty and second is that propagated from the determination of the lattice spacing. 





m K 


m N 


m A 


m s 


Win 


155.8 


553.7 


932.1(78.3)(14.4) 


1139.9(20.7)(17.6) 


1218.4(21. 5)(18.8) 


1393.3(6.7)(21.5) 


* 295.7 


593.5 


1093.1(18.9)(16.9) 


1253.8(14.1)(19.4) 


1314.8(15.4)(20.3) 


1447.7(10.0)(22.3) 


* 384.4 


581.4 


1159.7(15.4)(17.9) 


1274.1(9.1)(19.7) 


1316.5(10.4)(20.3) 


1408.3(7.0)(21.7) 


* 411.2 


635.0 


1214.7(11.5)(18.7) 


1350.4(7.8)(20.8) 


1400.2(8.5)(21.6) 


1503.1(6.5)(23.2) 


569.7 


713.2 


1411.1(12.2)(21.8) 


1503.8(9.8)(23.2) 


1531.2(11. 1)(23.6) 


1609.5(9.4)(24.8) 


701.4 


789.0 


1583.0(4.8)(24.4) 


1643.9(5.0)(25.4) 


1654.5(4.4)(25.5) 


1709.6(5.4)(26.4) 



LHPC [7] 

The LHPC collaboration calculates the light hadron spectrum in full QCD using a mixed 
action that exploits the lattice chiral symmetry provided by domain wall valence quarks 
(the D WF valence quark action) and ensembles of computationally economical improved 
staggered sea quark configurations (the so-called asqtad action). The lattice spacing is 
determined to be a = 0.12406(248) fm and the lattice volume is 20 3 x 64 . The range of 
pion masses simulated in this work extends from 758 MeV down to 293 MeV. 

Table 10. Masses of the pseudoscalar mesons and the octet baryons (in units of MeV) obtained by 
the LHPC collaboration (TABLE II, TABLE VI, TABLE VII of Ref. [7]). The uncertainties have 



the same 


origin as 


those in Table 9. 








m v 


m K 


m N 


m\ 


m s 


mn 


292.9 


585.6 


1098.9(8.0)(22.0) 


1240.5(4.8)(24.8) 


1321.6(6.4) (26.4) 


1412. 2(3. 2)(28.2) 


* 355.9 


602.9 


1157.8(6.4) (23.1) 


1280.2(4.8)(25.6) 


1350.2(4.8)(27.0) 


1432.9(3.2)(28.6) 


* 495.1 


645.4 


1288.2(6.4)(25.8) 


1369.3(4.8)(27.4) 


1409.1(6.4)(28.2) 


1469.5(4.8)(29.4) 


596.7 


685.6 


1394.8(6.4) (27.9) 


1440.9(8.0)(28.8) 


1463.1(9.5)(29.2) 


1504.5(8.0)(30.1) 


687.7 


728.1 


1502.9(11. 1)(30.0) 


1528.3(9.5)(30.6) 


1536.3(9.5)(30.7) 


1557.0(9.5)(31.1) 
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HSC [8] 



The HSC collaboration uses a Symanzik-improved action with tree-level tadpole-improved 
coefficients for the gauge sector and the anisotropic clover fermion action for the fermion 
sector. The lattice spacings are a s = 0.1227(8) fm and at = 0.003506(23) fm in spatial and 
temporal directions, respectively. The simulations are performed at four different lattice 
volumes 12 3 x 96, 16 3 x 96, 16 3 x 128, and 24 3 x 128. The simulated pion masses range 
from 383 MeV to 1565 MeV. For our purposes, we only need those data with m% < 0.5 
GeV 2 . The pseudoscalar meson masses and corresponding octet baryon masses are listed 
in Table 11. 



Table 11. Masses of the pseudoscalar mesons and the lowest-lying baryons (in units of MeV) 
obtained by the HSC collaboration (TABLE VI and TABLE VII of Ref. [8]). The uncertainties 
have the same origin as those in Table 9. 





m K 


m N 


m\ 


ms 


mn 


* 383.2 


543.6 


1147.5(10.7)(7.5) 


1243.1(8.4)(8.2) 


1287.0(8.4)(8.4) 


1347.8(6.8)(8.8) 


388.9 


545.9 


1164.9(22.5)(7.6) 


1226.8(16.9)(8.0) 


1288.7(16.9)(8.5) 


1345.0(11.3)(8.8) 


* 448.5 


580.8 


1238.1(16.9)(8.1) 


1328.1(11.3)(8.7) 


1361.9(16.9)(8.9) 


1412.5(10.7)(9.3) 


560.5 


646.6 


1361.9(22.5)(8.9) 


1440.6(16.9)(9.5) 


1457.5(22.5)(9.6) 


1496.9(16.9)(9.8) 



QCDSF-UKQCD [10] 

The QCDSF-UKQCD collaboration employs the particular clover action, which has a single 
iterated mild stout smearing, and the (tree-level) Symanzik improved gluon action, which 
contains the gluon action and the three-flavor Wilson-Dirac fermion action. The simulations 
are carried out at the lattice spacing of a ~ 0.075 — 0.078 fm, on 16 3 x 32, 24 3 x 48 and 
32 3 x 64 lattices. The resulting pion masses range from 229 MeV to 449 MeV. In the 
simulations, their kept the singlet quark mass fixed and tuned the quark masses to ensure 
that the kaon always has a mass less than the physical one. It should be noted that in 
Table 12 we did not tabulate the 16 3 x 32 and the 32 3 x 64 three-flavor simulation results, 
which have meson masses out of the range specified above. 

NPLQCD [11] 

The NPLQCD collaboration mainly studied finite-volume corrections on the octet baryon 
masses. Simulations are performed with rif = 2 + 1 flavor anisotropic clover Wilson action 
in four lattice volumes with spatial extent L ~ 2.0, 2.5, 3.0 and 3.9 fm. The anisotropic 
lattice spacing in the spatial direction is b s ~ 0.123 fm and bt = 6 s /3.5 in the time direction. 
The pion mass is fixed at ~ 390 MeV. 
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Table 12. Masses of the pseudoscalar mesons and the lowest-lying baryons (in units of MeV) 
obtained by the QCDSF-UKQCD collaboration (TABLE XX, TABLE XXII and TABLE XXIII of 



Ref. [10]). 


The uncertainties have the 


same origin as those 


in Table 9. 






m K 


m N 








228.9 


476.0 


999.2(13.6)(3.8) 


1144.2(15. 0)(4.3) 


1184.6(5. 8)(4.5) 


1266.6(9. 2)(4.8) 


97^ 9 


463.2 


1035.0(8.9)(3.9) 


1 1 97 fi"\ I A 
LIZ 1 .D^lZ.oJ^4.oJ 


1 1 79 (\S A P.M A A\ 


1 9/1 Q "i(Q 0^ I A 7\ 


* 320.4 


450.6 


1061.4(5.5)(4.0) 


1130.9(6.7)(4.3) 


1168.5(2.3)(4.4) 


1220.3(4.6)(4.6) 


324.2 


449.5 


1071.8(6.4)(4.0) 


1154.7(8.1)(4.4) 


1162. 8(2. 3)(4.4) 


1216.8(4.6)(4.6) 


350.7 


439.4 


1097.5(3.7) (4.1) 


1136.9(8.9)(4.3) 


1159.3(2.3) (4.4) 


1193.8(3.5)(4.5) 


383.0 


425.5 


1123.5(3.8)(4.2) 


1137.1(9.7)(4.3) 


1155.9(1.2)(4.4) 


1170.8(3.4)(4.4) 


411.5 


411.5 


1150.1(3.8)(4.3) 


1150.1(3.8)(4.3) 


1150.1(3.8)(4.3) 


1150.1(3.8)(4.3) 


448.8 


391.9 


1188.1(2.3)(4.5) 


1160.5(6.9)(4.4) 


1144.2(1.5) (4.3) 


1117.6(3.0)(4.2) 



Table 13. Masses of the pseudoscalar mesons and the lowest-lying baryons (in units of MeV) 
obtained by the NPLQCD collaboration (TABLE II of Ref. [11]). The uncertainties have the same 
origin as those in Table 9. 





m K 


m N 


ra A 


rriY, 


mn 


387.8 


544.4 


1182.1(5.4)(7.7) 


1263.3(5.1)(8.2) 


1286.6(4.3)(8.4) 


1361.5(4.1)(8.9) 


* 387.8 


544.4 


1164.0(3.2)(7.6) 


1252.0(2.6)(8.2) 


1280.5(3.0)(8.3) 


1356.4(2.6)(8.8) 


* 387.8 


544.4 


1151.6(2.5)(7.5) 


1242. 3(2. 6)(8.1) 


1282.7(2.2)(8.4) 


1349.3(2.1)(8.8) 


* 387.8 


544.4 


1151.3(2.6)(7.5) 


1241. 2(2. 2)(8.1) 


1279.0(2.8)(8.3) 


1349.2(2.0)(8.9) 



B Linear fits of M| = a + bM 2 

In this section, we determine the coefficients of the linear functions, M\ = a + bM 2 , by 
fitting the lattice data from the PACS-CS, LHPC, QCDSF-UKQCD and HSC collabora- 
tions. The fitting results are displayed in Fig. 6. It should be noted that the data point 
of the PACS-CS collaboration denoted by a red hollow circle is not included in the cor- 
responding fit because it has a lower strange-quark mass. The explicit expressions of the 
fitting functions are: 

• PACS-CS: 

Af£ = 0.291751 + 0.670652A/2. (B.l) 

• LHPC: 

M 2 K = 0.301239 + 0.479545M2 . (B.2) 

• QCDSF-UKQCD: 

M 2 K = 0.252658 - 0.489594M^ . (B.3) 
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Figure 6. (Color online) M% vs. M%. The dashed line is the fit to the PACS-CS data, hashed line 
the fit to the LHPC data, dotted line the fit to the QCDSF-UKQCD data, and dot-dashed line the 
fit to the HSC data. In the fit of the PACS-CS data, the ensemble with a lower strange-quark mass 
(the red hollow circle) is not included. 



HSC: 



M\ = 0.187873 + 0.734493M^. 



(B.4) 



- 23 - 



References 

[1] K. G. Wilson, Confinement of quarks, Phys. Rev. D 10 (1974) 2445. 

[2] C. Gattringcr and C. B. Lang, Quantum Chromodynamics on the Lattice - An Introductory 
Presentation. Springer, Heidelberg, 2010. 

[3] C. Alexandrou et al., Low-lying baryon spectrum with two dynamical twisted mass fermions, 
Phys. Rev. D 80 (2009) 114503, [arXiv : 0910 . 2419]. 

[4] S. Diirr et al., Ab-initio Determination of Light Hadron Masses, Science 322 (2009) 1224, 
[arXiv: 0906. 3599]. 

[5] S. Aoki et al. (PACS-CS Collaboration), 2+1 flavor lattice QCD toward the physical point, 
Phys. Rev. D 79 (2009) 034503, [arXiv : 0807 . 1661]. 

[6] S. Aoki et al. (PACS-CS Collaboration), Physical point simulation in 2+1 flavor lattice 
QCD, Phys. Rev. D 81 (2010) 074503, [arXiv : 0911 . 2561]. 

[7] A. Walker-Loud el al., Light hadron spectroscopy using domain wall valence quarks on an 
Asqtad sea, Phys. Rev. D 79 (2009) 054502, [arXiv: 0806. 4549]. 

[8] H.-W. Lin et al. (HSC Collaboration), First results from 2+1 dynamical quark flavors on an 
anisotropic lattice: light-hadron spectroscopy and setting the strange-quark mass, Phys. Rev. 
D 79 (2009) 034502, [arXiv : 0810 . 3588]. 

[9] W. Bictenholz et al., Tuning the strange quark mass in lattice simulations, Phys. Lett. B 690 
(2010) 436, [arXiv: 1003. 1114]. 

[10] W. Bictenholz et al. (QCDSF-UKQCD Collaboration), Flavor blindness and patterns of 
flavor symmetry breaking in lattice simulations of up, down, and strange quarks, Phys. Rev. 
D 84 (2011) 054509. 

[11] S. R. Beane et al. (NPLQCD Collaboration), High statistics analysis using anisotropic clover 
lattices: IV. Volume dependence of light hadron masses, Phys. Rev. D 84 (2011) 014507, 
[arXiv: 1104.4101]. 

[12] J. Beringer et al. (Particle Data Group), The Review of Particle Physics, Phys. Rev. D 86 
(2012) 010001. 

[13] Z. Fodor and C. Hoelbling, Light Hadron Masses from Lattice QCD, Rev. Mod. Phys. 84 
(2012) 449, [arXiv: 1203.4789]. 

[14] S. Weinberg, Phenomenological Lagrangians, Physica A 96 (1979) 327. 

[15] J. Gasser and H. Leutwyler, Chiral perturbation theory to one loop, Ann. Phys. (N. Y.) 158 
(1984) 142. 

[16] J. Gasser and H. Leutwyler, Chiral perturbation theory: Expansions in the mass of the 
strange quark, Nucl. Phys. B 250 (1985) 465. 

[17] J. Gasser, M. E. Sainio, and A. Svarc, Nucleons with chiral loops, Nucl. Phys. B 307 (1988) 
779. 

[18] H. Leutwyler, Principles of Chiral Perturbation Theory, Lecture (1994) [hep-pri/9406283]. 

[19] V. Bernard, N. Kaiser, and U.-G. Mcifincr, Chiral Dynamics in Nucleons and Nuclei, Int. J. 
Mod. Phys. E4 (1995) 193, [hep-ph/9501384]. 

[20] A. Pich, Chiral Perturbation Theory, Rep. Prog. Phys. 58 (1995) 57, [hep-ph/9502366]. 



- 24 - 



[21] G. Ecker, Chiral perturbation theory, Prog. Part. Nucl. Phys. 35 (1995) 84, 
[hep-ph/9501357]. 

A. Pich, Effective Field Theory, hep-ph/9806303. 

V. Bernard and U.-G. Meifiner, Chiral Perturbation Theory, Annu. Rev. Nucl. Part. Sci. 57 

(2007) 33, [hep-ph/0611231]. 

V. Bernard, Chiral perturbation theory and baryon properties, Prog. Part. Nucl. Phys. 60 

(2008) 82. 

S. Scherer and M. R. Schindler, A Primer for Chiral Perturbation Theory. Springer, 
Heidelberg, 2012. 

E. E. Jenkins and A. V. Manohar, Baryon chiral perturbation theory using a heavy fermion 
lagrangian, Phys. Lett. B 255 (1991) 558. 

T. Bccher and H. Lcutwyler, Baryon chiral perturbation theory in manifestly Lorentz 
invariant form, Eur. Phys. J. C 9 (1999) 643. 

J. Gegelia and G. Japaridzc, Matching Heavy Particle Approach to Relativistic Theory, Phys. 
Rev. D 60 (1999) 114038, [hep-ph/9908377]. 

T. Fuchs, J. Gegelia, G. Japaridze, and S. Scherer, Renormalization of relativistic baryon 
chiral perturbation theory and power counting, Phys. Rev. D 68 (2003) 056005, 
[hep-ph/0302117]. 

V. Bernard, T. R. Hcmmcrt, and U.-G. Meifiner, Cutoff schemes in chiral perturbation theory 
and the quark mass expansion of the nucleon mass, Nucl. Phys. A 732 (2004) 149, 
[hep-ph/0307115]. 

R. D. Young, D. B. Leinweber, and A. W. Thomas, Convergence of Chiral Effective Field 
Theory Chiral Effective, Prog. Part. Nucl. Phys. 50 (2003) 399. 

D. B. Leinweber, A. W. Thomas, and R. D. Young, Physical Nucleon Properties from Lattice 
QCD, Phys. Rev. Lett. 92 (2003) 24, [hep-lat/0302020]. 

R. D. Young and A. W. Thomas, Octet baryon masses and sigma terms from an SU(3) chiral 
extrapolation, Phys. Rev. D 81 (2010) 014503, [arXiv : 0901 . 3310]. 

A. Semke and M. F. M. Lutz, Baryon self energies in the chiral loop expansion, Nucl. Phys. 
A 778 (2006) 153, [nucl-th/0511061]. 

E. E. Jenkins, Baryon masses in chiral perturbation theory, Nucl. Phys. B 368 (1992) 190. 

V. Bernard, N. Kaiser, and U.-G. Meifiner, Critical Analysis of Baryon Masses and 
Sigma-Terms in Heavy Baryon Chiral Perturbation Theory, Z. Physik C 60 (1993) 111, 
[hep-ph/9303311]. 

M. K. Banerjee and J. Milana, Baryon mass splittings in chiral perturbation theory, Phys. 
Rev. D 52 (1995) 11. 

B. Borasoy and U.-G. Meifiner, Chiral expansion of baryon masses and a-terms, Ann. Phys. 
(NY.) 254 (1996) 192, [hep-ph/9607432]. 

A. Walker-Loud, Octet Baryon Masses in Partially Quenched Chiral Perturbation Theory, 
Nucl. Phys. A 747 (2004) 476, [hep-lat/0405007]. 

P. J. Ellis and K. Torikoshi, Baryon Masses in Chiral Perturbation Theory with Infrared 
Regularization, Phys. Rev. C61 (1999) 015205, [nucl-th/9904017]. 



- 25 - 



M. Frink and U.-G. Meifiner, Chiral extrapolations of baryon masses for unquenched 
three-flavor lattice simulations, J HEP 07 (2004) 028, [hep-lat/0404018]. 

M. Frink, U.-G. Meifiner, and I. Scheller, Baryon masses, chiral extrapolations, and all that, 
Eur. Phys. J. A. 24 (2005) 395, [hep-lat/0501024]. 

B. C. Lchnhart, J. Gegelia, and S. Scherer, Baryon masses and nucleon sigma terms in 
manifestly Lorentz-invariant baryon chiral perturbation theory, J. Phys. G: Nucl. Part. Phys. 
31 (2004) 89, [hep-ph/0412092]. 

J. Martin Camalich, L. S. Geng, and M. J. Vicente Vacas, Lowest-lying baryon masses in 
covariant SU(3)-flavor chiral perturbation theory, Phys. Rev. D 82 (2010) 074504, 
[arXiv: 1003. 1929]. 

A. Semke and M. F. M. Lutz, On the possibility of a discontinuous quark-mass dependence of 
baryon octet and decuplet masses, Nucl. Phys. A 789 (2007) 251. 

A. Semke and M. F. M. Lutz, Quark-mass dependence of the baryon ground-state masses, 
Phys. Rev. D 85 (2012) 034001. 

P. C. Bruns, L. Greil, and A. Schafer, Chiral extrapolation of baryon mass ratios, 
arXiv: 1209.0980. 

M. F. M. Lutz and A. Semke, On the consistency of recent QCD lattice data of the baryon 
ground-state masses, arXiv : 1209 . 2791. 

K. I. Ishikawa et al. (PACS-CS Collaboration), SU(2) and SU(3) chiral perturbation theory 
analyses on baryon masses in 2+1 flavor lattice QCD, Phys. Rev. D 80 (2009) 054502, 
[arXiv: 0905. 0962]. 

J. M. Alarcon, L. S. Geng, J. Martin Camalich, and J. A. Oiler, On the strangeness content 
of the nucleon, arXiv : 1209 . 2870. 

J. A. Oiler, M. Verbeni, and J. Prades, Meson-baryon effective chiral Lagrangians to 0(q 3 ), 
JEEP 09 (2006) 079. 

L. S. Geng, J. Martin Camalich, L. Alvarez-Ruso, and M. J. Vicente Vacas, Leading 
SU (3) -Breaking Corrections to the Baryon Magnetic Moments in Chiral Perturbation 
Theory, Phys. Rev. Lett. 101 (2008) 222002. 

L. S. Geng, M. Altenbuchinger, and W. Wcisc, Light quark mass dependence of the D and 
decay constants, Phys. Lett. B 696 (2011) 390. 

L. S. Geng, J. Martin Camalich, L. Alvarez-Ruso, and M. J. V. Vacas, Lowest lying spin- 1/2 
and spin-3/2 baryon magnetic moments in chiral perturbation theory*, Chinese Phys. C 34 
(2010) 1307, [arXiv: 1001.0465]. 

J. Bijnens and I. Jemos, A new global fit of the L\ at next-to-next-to-leading order in Chiral 
Perturbation Theory, Nucl. Phys. B 854 (2012) 631-665, [arXiv : 1103 . 5945]. 

L. S. Geng, X.-L. Ren, J. Martin Camalich, and W. Weise, Finite-volume effects on 
octet-baryon masses in covariant baryon chiral perturbation theory, Phys. Rev. D 84 (2011) 
074024. 

G. Amoros, J. Bijnens, and P. Talavera, QCD isospin breaking in meson masses, decay 
constants and quark mass ratios, Nucl. Phys. B 602 (2001) 87. 

L. S. Geng, X.-L. Ren, J. Martin Camalich, and J. Meng, Finite-volume effects on 
octet-baryon masses in covariant baryon chiral perturbation theory revisited, in preparation. 



- 26 - 



[59] M. Procura, T. R. Hcmmcrt, and W. Weise, Nucleon mass, sigma term and lattice QCD, 
Phys. Rev. D 69 (2003) 034505, [hep-lat/0309020]. 

[60] M. Procura, B. U. Musch, T. Wollenweber, T. R. Hemmert, and W. Weise, Nucleon mass: 
from lattice QCD to the chiral limit, Phys. Rev. D 73 (2006) 114510, [hep-lat/0603001]. 

[61] J. Giedt, A. W. Thomas, and R. D. Young, Dark matter, the CMSSM and lattice QCD, 
Phys. Rev. Lett. 103 (2009) 201802, [arXiv : 0907 . 4177]. 

[62] S. Diirr et al. (BMW Collaboartion) , Sigma term and strangeness content of octet baryons, 
Phys. Rev. D 85 (2012) 014509. 

[63] P. E. Shanahan, A. W. Thomas, and R. D. Young, Sigma terms from an SU(3) chiral 
extrapolation, arXiv : 1205 . 5365. 

[64] A. Semke and M. F. M. Lutz, Strangeness in the baryon ground states, arXiv : 1202 . 3556. 

[65] R. Horsley et al. (QCDSF Collaboration), Hyperon sigma terms for 2+1 quark flavors, Phys. 
Rev. D 85 (2012) 034506. 

[66] J. Gasser, H. Leutwyler, and M. E. Sainio, Sigma term update, Phys. Lett. B 253 (1991) 252. 

[67] J. M. Alarcon, J. Martin Camalich, and J. A. Oiler, The chiral representation of the ttN 
scattering amplitude and the pion-nucleon sigma term, Phys. Rev. D85 (2012) 051503, 
[arXiv: 1110.3797]. 

[68] G. Bali et al. (QCDSF Collaboration), Strange and light quark contributions to the nucleon 
mass from lattice QCD, Phys. Rev. D 85 (2012) 054502. 

[69] G. Bali et al. (QCDSF Collabortion), Nucleon mass and sigma term from lattice QCD with 
two light fermion flavors, arXiv : 1206 . 7034. 

[70] H. Ohki et al. (JLQCD Collaboration), Nucleon sigma term and strange quark content in 
2+1-flavor QCD with dynamical overlap fermions, Proceedings of Science LAT2009 (2009) 
124, [arXiv: 0910. 3271]. 



-27- 



